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Abstract 

We study the dynamics of a classical particle moving in a punctured 
plane under the influence of a strong homogeneous magnetic field, an 
electrical background, and driven by a time-dependent singular flux tube 
through the hole. 

We exhibit a striking classical (de)localization effect: in the far past 
the trajectories are spirals around a bound center; the particle moves 
inward towards the flux tube loosing kinetic energy. After hitting the 
puncture it becomes "conducting": the motion is a cycloid around a 
center whose drift is outgoing, orthogonal to the electric field, diffusive, 
and without energy loss. 

PACS numbers: 45.50. Pk Particle orbits classical mechanics, 45.50.-j Dy- 
namics and kinematics of a particle and a system of particles, 73.43.-f 
Quantum Hall effects, 73.50.Gr Charge carriers: generation, recombina- 
tion, lifetime, trapping, mean free paths 



1 Introduction 

The motivation to study the dynamics of this classical system is to sharpen 
our intuition on its quantum counterpart which is, following Laughlin's and 
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Halperin's [12] proposals, widely used for an explanation of the Integer Quantum 
Hall effect. Of special interest is how the topology influences on the dynamics. 
In the mathematical physics literature Bellissard et al. (5] and Avron, Seller, 
Simon ^ used an adiabatic limit of the model to introduce indices. The 
indices explain the quantization of charge transport observed in the experiments 
[13] . See [3 Uni IHl ISl HI] for recent developments. We discussed the adiabatics 
of the quantum system in [2], its quantum and semiclassical dynamics will be 
treated elsewhere. The dynamics of the classical system without magnetic field 
were discussed in yj. 

We state the model and our main results: 

Consider a classical point particle of mass m > and charge e > moving 
in the punctured plane \ (0). Suppose that a magnetic flux line with time 
varying strength $ pierces the origin and further the presence of a homogeneous 
magnetic field of strength B > orthogonal to the plane and an interior electric 
field with smooth bounded potential V. 

The equations of motions are Hamiltonian. For a point 
{(1,P) = ((^1,92), (Pi,P2)) in phase space 

P = \ (0) X 
the time dependent Hamiltonian is : 



where := (— q'2,9i)- We suppose that 

$ : R ^ M and V : M X ^ R are smooth functions. 
The electric field is —dtA — dqV , the force on the particle with velocity q: 

e {q A rotiA) - dtA - d^V) = -e (^Bq^ - MlJ^ + d,V^ 

Remark that the part of the electric field induced by the flux has circulation 
but vanishing rotation, and is long range with an 1/r singularity at the 
origin, we call it the circular parts. V is smooth on the entire plane so that the 
circulation of the corresponding field is zero. This is the topology essential for 
the dynamics. 
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Recall that when only the constant magnetic field is present, the particle 
follows the Landau orbits; these are circles around a fixed center with frequency 
^ whose squared radius is proportional to the energy. 

Our result for the case $ ~ t, i? large, V such that the torque q A dgV is 
small is qualitatively: 

- the motion in configuration space is approximately rotation with radius 
proportional to the square root of the (time-dependent) energy around a 
drifting center. 

- for large enough negative times the center is trapped by the flux line and the 
energy is linearly decreasing with time, so the particle is spiraling inwards 

- from the hitting time on (i.e. the time when the Landau orbit "hits" the 
singularity) the center starts to drift away from the flux line, the energy 
remains asymptotically constant in the future. The drift is diffusive. The 
situation is described by Fig.[T] showing a typical orbit in q-space. 

Remark that the corresponding analysis remains true if the sign of B is changed. 

In this case we may state our observation as: Hall conducting states are eventually 

trapped by the flux line and trapped states are energy conducting. 

Here "hall conducting" means that the center follows the lines of the potential 

diffusively. 

We shall discuss the corresponding quantum behavior elsewhere. 

In the first section of this paper we state some general remarks on the model 
and discuss the problem for frozen values of the flux. Next we define appropriate 
action angle coordinates and use an averaging (adiabatic) method to approximate 
the dynamics near the hitting time between the particle and the flux line. In 
the last section we discuss the asymptotic behavior of the solution of the full 
equations of motion. 

Let us remark that our method includes (for the two dimensional case) a 
simple proof for the guiding center approximation widely used in plasma physics. 



2 Dynamics of the frozen system 



Denote 
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Figure 1: Typical trajectory of the Hamiltonian ^ — y^Q^ — + sdqVjj with 
V chosen to be V{x,y) = l/3(sina; + siny) 
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We use the scaling (t, q, p) i— ^ {cut, q/ A, p\) and "absorb " V into the time depen- 
dent vectorpotential. The scaled variables are called (s,q,p). The Hamiltonian 
under consideration then reads 

H{s;p, q) := ^ (p - a(s; q)f ; a(s; q) := Qg^ + a£;(s; g)^ 

where a£;(s) : \ (0) ^ is smoothly time dependent with rot{aE){s) = 0. 
aE{s) and the electric field E{s) : \ (0) ^ are defined by: 

-S,M.):^EW:^i(§i(l)^-A(W(5.A,)) (1) 

We discuss first the solution of the equation of motions for a frozen time 
(T G M. As dsaE{(T] q) = 0, the solution of the frozen equations generated by the 
Hamiltonian H(a) goes along the lines of the classical Landau problem (which 
means: the case $ = 0; V{q) = 0) 

For a G M define 

1. the velocity field: v{a) : P M^, v{a; q,p) := p — a{a; q); 

2. the center: c(cr) : P R^, c{a; q,p) := q — v^{a; q,p); 

3. the angular momentum: L : P ^ M, L{q,p) := q Ap. 
Denote the Poisson bracket: {f,g} = dgfdpg — dpfdqg. 

We list some useful formulas: 

Proposition 2.1 The following identities hold as functions on phase space P for 

all a G M; 

1. {Vi,V2} = l, {Ci,C2} = -l, |c,^j|=C-L, {ci,Vj} = 0; 

2. H= {v, H} = -v^, {c, H} = 0; 
3. 

1 1 

-c^ = -v^ + L- qAaE = H + L- qAaE] (2) 
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4. the frozen flow {q{(T; s),p{a; s)) defined by 
dsq{a; s) = dpH{a), dsp{a; s) = -dqH{a), 
(g(a;0),p(a;0)) = {q,p) is : 

q{a; s) = c{a) + cos{s)v'^{a) + sm{s)v{a) 

p{a; = ^ + cos(s)f (a) - sin(s)f "^(a)) + a£;(a; q{a; s)) 

Proof (1),(2),(3): {vi,V2} = {pi - ai{a, q) , p2 - a2{a, q)} = rot{a{a)) = 
1, {gi, Vj} = 5ij. H = ^v"^ so {g, H} = v, {v, H} = -v^. = q^ + v'^ + 2qAv; 
on the other hand L = q /\v + \q^ + g A a£;(cr; q). 

(4): The force is —q-^ independently of a, Newton's equation g = — g-*- is readily 
verified. On the other hand: p = v + a = a + c^ — q-^ = c^ — \q-^ + a£;((T; g). 
So p{s) follows from g(s) □. 

Remarks 2.2 1. Since the energy H{a) = (cr)^ is conserved under the 
frozen flow, the projections of the trajectories to q-space are circles around 
c(cr) with radius ^/2H{a). An orbit encircles the origin (has non-trivial 
homotopy) in \ (0) // and only if 

< 2H L - q A aE{cT] q) < 0; 

2. the flow is, strictly speaking, not complete as for L — q A aEicr; q) = 
the particle reaches the origin in q-space (and infinity in p-space) in finite 
time; the energy remains, however, finite. This is a mathematical subtlety 
which can be handled. 



3 Action angle coordinates 

in order to discuss the full dynamics for large B we introduce action angle co- 
ordinates. The frozen dynamics as discussed in Proposition 12. II suggests to take 
as coordinates the angles and absolute values of c and v^, i.e. with the 

notation: e{9) := {cos9,sm6) : 



q = c + V 



+ V 



--: |c|e(v3i) + \v\e{-ip2) 



1 1 

p = - (c^ + f) + asia; q) = - [\c\e^{^i) - |t'|e^(-^2)) + aE^a; g) 



Motivated by this we define for a G M 

q{a;Lp,I) := ^/2he{ipi) + ^/2I^e{-ip2) 

p{a; if, I) := - (^^/2he^ (ifi) - ^/2he^{-ip2)^ + asia; q{a; ip, I)) 

and, denoting by C the nullset 1)](pi + Lp2 = tt, h = h} where g(cr; (p, /) = 
0, by V the nullset f ^ = or = 0}. Thus for each frozen time a G M 

the transformation to action angle coordinates T(cr) is defined by 

T(a) : 5^ X 51 X {(/i, Is); Ji > 0, /2 > 0} \ C ^ P \ P 
T(o-;v?,/) = T(o-;<^i,v92,/i,/2) := (g(o-; v?, /)) 



We have 

Lemma 3.1 1. T{a) is a canonical diffeomorphism 
2. T^^{a) is determined by 



c^(a) 1 / / I 



J2(a) = = - + g)^ ^ 



c 



v/2(iJ(f7) + L-gAas((T; g)) 



e{-ip2ia)) = j- 



lq + P^ + a^{a]q) 



Proof: These identities follow immediately from Proposition 12.11 

{/i, I2} = 0, {e(<^i), e{ip2)} = 0, {/i, e(<^2)} = = {h, e(<^i)}, 

1 c-^ 
{e((^i), /i} = --{c, — } = — = e^{ipi). 

On the other hand, {e{(pi),Ii} = e^{(pi){(pi, h}, so {tpi, Ii} = 1. Similarly: 

{(^2,/2} = l. □ 
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We now investigate the full equations of motion, i.e. those for time-dependent 
flux, in these action angle coordinates. As rot(E) = there exists a (possibly 
multi-valued) function which we denote by m = m{s;q) such that 

dqm{s) = E{s) = -dsaE{s). 

Then T(s) is generated by m: 

dsT{s; LP, I) = (0, dsaEiqis; /)) = {dpm, -dgm) o T{s; (p, I). 

Denote by U{s) : P ^ P the hamiltonian flow of H{s) defined by U{s) : = 
{q{s),p{s)) 

q{s) = dpH, p{s) = -d,H, (g(0),p(0)) = {q,p), 
then for the flow U{s) = {(p{s), I{s)) in action angle coordinates defined by 

T{s)oU{s) = U{s)oT{s = 0) 

it holds: 

= djK o U{s), i{s) = -d^K o ((/.(O), J(0)) = J), 

where the Hamiltonian in action angle coordinates, K = HoT — moT, is 

K{s] (p, I) = I2- m(s; q{s; (/?, /)) 
and the equations of motion are (with the notation (•, ■) for the scalar product) 

= djK=(^^^^- {Eis, qis; ^, I)), djq) (3) 
i{s) = -d^K = {E{s; q{s; J)), d^q) (4) 

Remark 3.2 Another way to derive these equations is to start from Newton's 
equation 

q = -q^ + E{s;q). 
From the very definition of c and v one gets: 

c=-E^{c + v^) V = -v^ + E{c + v^) 

which in action angle coordinates gives (0), (0)- 
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4 Averaged dynamics 

We apply averaging with respect to the fast angle (/?2 to the system (0), (jlj) (see 
(TSlI^). The singularity problem can be overcome by a regularization technique 
(see [16j). The result is that the solutions of the equations are at a distance of 
order 1/5 to the solution of the averaged equations over times of order B. 

Remark that at this place we are mainly interested in the (de)localization 
effect so we did not make use of more involved adiabatic or KAM methods in 
order to go to longer or even infinite time scales. 

We detail this for the case 



In particular for a function / defined on the plane thus depending only on the 
variable q we denote 



^(t) = $ot, V time independent, 



i.e., a flux $o per unit time is added ad-eternam. 

Denote the average of a function / on the phase space by 





The field dU is 




Define 



f 



e$o 



271UJ 



and choose m and thus K: 



f arg{q) 



-V{\q) 
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Making use of the identities 
the system (jH) reads 



\ „ sin((y9i + (P2) ( y h \ e 



^ ^ V 1 / 2(/i + /2 + 2VJ^cos(^i + (^2)) 1 -./IT / ' ^ 




'^'^ = '2(/, + /, + 2'^cos(v.. + ^.)) ( 1 ) + K -1) - 5 
The averaged quantities are readily calculated: using 

1 \ 1 /sin {ifi + (^92) 



one finds for the averaged vectorfield 



>/2) \ e / d^,Va,{ipi,X^I) 



UJ 







where we used the binary function x'- x{True) := 1, x(-^o/se) := 0. 

Remark 4.1 Remark that the averaged vectorfield is the hamiltonian vectorfield 
derived from the from the "averaged" Hamiltonian Kav- Indeed, using the split- 
ting of arg{q), which is a multi-valued function defined on the covering space of 
\ (0), into a linear and oscillating part 



arg(g(</?,/)) 



^1 + arg(^(l,0) + y^e(-</?i -^2)) if h > h 
-V^2 + arg((l,0) + y/|^e(^i + ^2)) if h > h 
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and: 



2n 



arg((l,0) + ae(s))cis = for < a < 1; 



One finds that for 



one has d^Kav = {d^K)av, djKav = idjK)av 

The result on the dynamics now is: 

Theorem 4.2 Denote by J = (Ji, J2), ip = (^1,^2) the solution of the av- 
eraged equations (|3p 

ij{s) = djK,,{ij{s),J{s)), J(0) = (J0,J°) 

and by I = (/i, I2), (p = (</?i, (P2) the solution of the full equations (0) 

= djK{v{s),I{s)), /(0) = (/?,/°) 
lis) = -d^Ki^is),Iis)), ifiO) = {iflifl) 

then it holds 

1. LetV = 0, denote A J = J° - Ji then: 



J{s) = mm{J^, J2} + (fs - ^J) 



X(fs> AJ) 
-x{is<AJ) 



2. For any V and any si, S2 G M 

1^2(^2)- J2(si)|=f 



x{Ji{u) < J2{u))du 
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3. Let V be such that the torque of the corresponding field satisfies for a 

ce[0,l): 

\qAd,V\<% 
Ztt 

then for any initial condition it holds: 



Ii — I2 is strictly increasing, furthermore 

f (1 - c) < /i(s) - /2(s) < f (1 + c) (Vs e M). 
4. In particular if q A d(jV = it holds for all s e R: 

I,{s) - his) = i{s - So) (6) 

where sq 's the unique "hitting time" defined by this equation. 

Proof Using that for V = it holds Ji(s) - J2{s) - is = AJ the first 
assertion follows by inspection. The second assertion follows from integration of 
1^. Finally we have from l^: 

h-i2 = dsih-l2) = qAE = -(^-iqA d,V)i\q) 

from which the last assertion follows. □ 

Remarks 4.3 1. The first equation explains the qualitative behavior of the 
solution exhibited in Fig. d- Ji is linear in time in the future and is constant 
in the past. 

2. Loosely speaking the second assertion of the theorem means that, on the 
average, one has 

I energychange] = \ fluxchange through the orbit during stay time\ 

where the stay time means the time where the "orbit surrounds the origin". 
This should be like this as the the change in energy equals the work of the 
electric field along the orbit: 



H{s;q{s)) - H{so;q{so)) = / {aE{s),ds). 

J sn 
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3. The last assertion says that the orbit presented in Fig.\^in the introduction 
is generic, i.e.: inward spi rating motion with fixed center followed by the 
usual Hall cycloids with the center following the lines of the potential. We 
argue that our condition on the potential is far from optimal and that for 
large enough magnetic field the situation described in this paper is generic 
for V smooth and bounded with bounded derivative. This needs further 
investigation. 



5 Large time asymptotics, potential free case 

For the case = <l>o^. V" = we can determine the large time asymptotics of 
the solution. We keep the notation f := Observe also that 

K = K{ip, I) = h- arg(v^e((/.i) + v^e(-(/.2)) 
is an integral of motion. 

Theorem 5.1 Denote by I = (/i,J2), = {ipi,ip2) the solution of the full 
equations of motion (0), (0) 

0(.) = djK{^{s),I{s)), J(0) = (J?,J°) 

i{s) = -d^K{^{s),I{s)), m = i^i,^l) 

then the following asymptotic behavior holds: 
in the future, s — » oo 

The following limits exist and define the constants oq > 0, bo: 

lim J2(s) =: , lim (v9i(s) + v?2(s) - s) =: 6o, \im. {his) - iipi{s)) = K, 

s— »oo 4l s— »oo s— >oo 

the asymptotics are 

his) : 



h{s) = h[s) + f{s-so) 

4f2 f \^_53/2 

/ \ , al K i , , , 1 
</?2(s) = s + bo--^ + -r cos(s + 6o)^ 



4f2 \ 1 

-l + 2cos(2(s + 6o))-— sin(2(s + 6o)) - + 3/2 
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with So defined as in 
in tlie past, s ^ —oo 

Ttie following limits exist and define the constants ao > 0, bo: 

lim Ji(s)=:-J, lim {ipi{s) + (p2is) - s) =: bo, lim {l2{s) + i(p2is)) 

s— >— oo 41 s-^— oo s— »— oo 

the asymptotics are 



4f V 2 ^ ' ^7 ^ 4 V 2f ' ' s 

his) = h{s)-i{s-so) 

/ \ 7 K f , ~ . 1 

V5i(s) = So + Oo + - + — cos(s + 6o) 



4f2 f ao ' 

1 / ~ 4f2 ~ \ 1 

1-2 cos(2(s + 6o)) - — sin(2(s + bo)))- + 



s \\s\y^ 



^2{s) = s-so--^ + -r--r + 



4f2 ' f 4s ' \^|5|3/2 



Proof We give an outline of the main steps of the proof for the case t ^ oo. 
Some particular computations in the proof turned out to be quite tedious and 
thus computer algebra systems were employed to facilitate them. 

Suppose t > 

Step 1 

From d^I) we know — his) = f(s — sq). So the equations of motion 
only involve J := Ji + h and ip := ipi + and transform to 

i^ssinib ■ f^s 

ih = 1 + — , , J 



V.P - Ps^J + VJ^ - f COS i;) J + VJ^-Ps^ cosV 

Step 2 

Do a second transformation 



Xi := V'P — f^s^ COS'?/', X2 := \J .P — f^s^ sin-?/', 
the J,?/' equations transform to 

x\ h X2 = -F(s,xi,X2), ±2 - a:i = 0, 
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with 

F{s, Xi,X2) := f 



S ^X\ + {X2 - f)2 + f2s2 + ■ 

The corresponding homogeneous system is equivalent to 

xi / 1 \ 

x\ I-1H — ^ a^i = or si/ + ?/ + s?/ = 

s \ ) 

with y defined by x\ = sy. The latter is Bessel's equation of order so one has 
two independent solutions of the homogeneous system: 

xiis) \_( sMs) \ ( x,{s) \_( sYo{s) 

X2{S) ) \ SJ,{S) ) ^""'^ \ X2{S) ) V ^^l(^) 

with the Bessel functions Jm (^m) of the first (second) kind. 
Step 3 

Transform the x-differential equation to the integral equation 

xi{s) = cisJo{s) + C2sYo{s) 

ITS r°° 

- — j {Y,{s)Ut) - Us)Y,{r))F{T,x,{T),X2{r))dT 
X2{s) = CisJi{s) + C2sYi{s) 

-y y (n(s)Ji(r)- Ji(s)Fi(r))F(r,Xi(r),X2(r))dr 

where the numbers Ci,C2 involve the initial conditions. 

The equation is of the form x = K,{x), the solution is constructed as 
the limit of the sequence Xn+i = /C(x„) starting from xq = 0. To ver- 
ify the convergence one can apply yet another substitution x{s) = y{s)/y/s, 
G{s,y) = s~^/'^F{s, s~^^'^y). Consequently the integral equation takes the form 

/■oo 

y{s)=yo{s)- I J^{s,T)G{T,yi{T),y2{T))dT 

J s 



where 



yOj{s) = Cl^/sJj^l{s) + C2y/sYj^i{s), j = 1, 2, 

J^j{s,t) = ^v^(F,„i(s)Ji(r)- J,_i(s)Fi(r)), j = 1,2. 
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Considering the new integral equation in the Banach space L°°([ s*, cx)[) ® 
one can show that the iteration process is indeed contracting provided > 1 is 
sufficiently large. It is then straightforward to derive from the integral equation 
the asymptotic expansion of the solution x{s). One finds that 

/ 5 \ 1 /I 

= aoe{t + bo)V7s+ ( ^ e(t + 6o) - gOo e^(t + 6o) 1 + 

Step 4 

Transforming back first to the J, then to Ji, /2, v^i, </?2 variables gives the 
claimed asymptotic expansion. □ 

The asymptotic formulae for the actions and the angles imply the following 
asymptotic behavior of the solutions and the energy thus defining the transport 
coefficients: 

Denote 

the energy in the original coordinates q,p, and g^c = q/^ the scaled coordinate. 
Rescaling then gives 

Hit) = uH{ujt) = uhiujt), q{t) = Xqsc{ujt), q^c = ^/2he{Lpi) + ^/2I^e{-ip2). 

This leads to the following limits valid for any fixed initial condition and any 

B >0,%>0: 



qit) Jal K 
— -;=r — >f^nn \ / - — :r e — -r — 



V 2nB \AP f 

e[—ujt) 



2 



4f 
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